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I. INTRODUCTION 



Multipartite entanglement of three trapped ions in a cavity and W- State generation 
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A scheme to generate three qubit maximally entangled W-states, using three trapped ions in- 
teracting with red sideband tuned single mode field of a high finesse cavity, is proposed. For the 
cavity field initially prepared in a number state, the probability of generating three ion W-state is 
OO \ calculated. By using the ion-cavity coupling strengths achieved in experimental realizations, the 

>P . interaction time needed for W-state generation is found to be of the order of 10 /i sec. It is found 

that for a fixed number of photons in the cavity the nature of entanglement of ionic internal states 
can be manipulated by appropriate choice of initial state phonon number. 

The ionic qubits in W-like state are found to be entangled to cavity photons. Analytical expres- 
sions for global negativity and partial ii'— way negativities {K = 2 to 4) are obtained to study the 
evolution of entanglement distribution as a function of interaction parameter. Reversible entangle- 
C^ , ment exchange between different entanglement modes is observed. For specific values of interaction 

CN . parameter, the three ions and photon-phonon system are found to have four partite entanglement, 

generated by 2— way and 3— way correlations. 

-)— > 

^. 

^ O"^. Interest in manipulation of multipartite quantum states is motivated by the possible use in quantum communication 

[l| , quantum dense coding [2| and quantum teleportation Sj . Advances in experimental techniques of trapping, cooling, 
^sl • and manipulation of internal states of ions through interaction with external lasers [^, |5|, |6| have made the generation 
J> ' of maximally entangled states Q a reality. While single atoms and ions, with long-lived internal states, are suitable for 
^^ , storing quantum information, the photonic qubits serve as fast and reliable carriers to transport quantum information 
^^ ' over long distances. Cirac et al. [8| considered the possibility of using photonic channels to interconnect high finesse 
cavities with single trapped ions inside. It was shown in ref. Q, that the composite quantum system with a single 
trapped two-level ion interacting with a quantized light field in single-mode cavity evolves into maximally entangled 
\l ' three-particle Greenberger-Horne-Zeilinger state. A different scheme to generate three qubit maximally entangled 
^^ , GHZ state, using a trapped ion interacting with a resonant external laser and sideband tuned single mode of a 
cavity field has been proposed in ref. \\ff\ and the dccoherence process of this three qubit system examined ll]| . 
Experimentally successful efforts in this direction 12, 13, 14,], have demonstrated coherent coupling of electronic and 
motional states of a single trapped Ca"*" ion to single field mode of a high finesse cavity. In the experiment of ref. 
[12| . a Ca"*" ion in a RF-Paul trap was placed with high precision at an arbitrary position in the standing wave field 
«2j ■ of a near-confocal resonator for several hours of interaction time. Since the size of the cavity is large compared to 
C^ ■ the wavelength of standing wave, a large number of nodes and antinodes of the standing wave are formed inside the 
cavity. The electromagnetically induced transparency (EIT) method used in the work of Morigi et. al. [l^l and Roos 
et al. |16|, could make possible the initialization of motional states of ion strings in a cavity. These experiments open 
the possibility that more than one trapped ions may be cooled down to their lowest vibrational states and interact 
with the quantized cavity field. A string of ions in the same trap with the center of mass in a definite motional mode 
may, likewise, be prepared inside a cavity. In this article, we propose a scheme in which three two-level trapped 
ions are cooled inside a high finesse cavity to generate multipartite entangled states through coherent coupling of 
their electronic and vibrational states to the standing wave cavity field. For the cavity field initially prepared in a 
number state, we obtain conditions on the interaction parameters to generate states in which three ions in W like 
state are entangled to cavity photons. Tessier et al. [ij] have investigated the multipartite entanglement of two 
two-level atoms, generated by resonant coupling to a single mode of the electromagnetic field. The resulting physical 
system corresponds to the two-atom Tavis-Cummings model [l8[. Our system is similar to that of ref. 17] in that the 
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entanglement between the ions arises not through direct interaction between the ions, but through their couphng to 
the phonon-photon system. We have not considered the effects of spontaneous decay, decoherence of motional states 
and cavity losses in this article. The object of studying a closed system, without decoherence effects, is to have a 
clear picture of dynamics of quantum correlations. The entanglement of mixed multipartite states resulting due to 
decoherence is not so well understood. The memory qubits (ionic internal states) are known to have a coherence time 
of the order of seconds and even minutes (magnetically insensitive transitions). The success of experiments with a 
single ion indicates, that the motional decoherence caused by trap fluctuations or environment noise is negligible on 
the time-scale of state generation and photon qubit coherence of the order of milisecs. The decoherence caused by 
the implementation of the scheme itself needs to be studied theoretically for the three ion case though such studies 
for single ion case have been done in refs. [lll[ig|. 

Multipartite entanglement is a useful resource for implementation of quantum information processing protocols 
[2(j |. To retrieve information about the entanglement of a group oi K {2 < K < 4) subsystems of the composite 
system, we use the global negativity and the partial K—way negativities. For analyzing the entanglement dynamics, 
we have considered the vibrational phonons and cavity photons to constitute a single quantum system. Negativity 
[2l|, l22j 1231 , based on Peres Horodecki [2J, [25| criterion, has been shown to be an entanglement monotone 26] . The 
K—way negativity refers to the negativity of a partial transpose constructed by imposing specific constraints during 
transposition. The coherences of a multipartite composite system having N subsystems can be quantified by K—way 
negativities [23| {2 < K < N). Partial K—way negativity is the contribution of a specific K—way partial transpose to 
global negativity For canonical states, the partial K—way negativities measure the genuine K— partite entanglement 
of the system [2^[23]- Genuine isT— partite entanglement refers to iiT— partite entanglement due to correlations simlar 
to those present in a GHZ-like state of K subsystems. Loss of a single subsystem destroys this type of entanglement 
completely, leaving no residual entanglement. We have obtained analytical expressions for partial K—way negativities 
with respect to photon-phonon state and internal states of ions. The advantage of using partial K—way negativities 
lies in the fact that the entanglement between parts of a composite system is obtained from the full state operator, 
without state reduction. With the cavity prepared initially in a Fock state, the entanglement of ionic internal states 
depends strongly on the choice of phonon number at t = 0. Extra control on qubit state manipulation, gained by 
coupling the vibrational modes to the cavity field |13l . Il4| , is an advantage for successful implementation of quantum 
gates. 

The vector space and system Hamiltonian are discussed in section II. Analytical expressions for the state of the 
system at current time, with the three ions prepared in their respective (i) ground states, and (ii) excited states, are 
presented in section III. The role of initial state center of mass quanta in state manipulation and the time evolution 
of W-state generation probabilities is also discussed in section III. In section IV, the entanglement dynamics of ions, 
phonons and photons is investigated. The Global and partial K—way negativities, used to study the entanglement 
distribution in the composite system, are defined and calculated analytically as well as numerically for special initial 
state preparations of the system. We discuss, briefly, the information gained from the dynamics of global and partial 
K—way negativities in section V, followed by the conclusions in section VI. 

II. THE MODEL 

Consider three two level trapped cold ions vibrating with trap frequency v inside a high finesse cavity. The frequency 
of standing wave cavity field is Uc and ions are well separated from each other so that no dipole interaction takes 
place between the ions. The free Hamiltonian of the system composed of internal states of ions, the vibrational state 
of center of mass and the state of standing wave cavity field is given by 

Ho = hv U^a + i) + hwM +^J2a^^\ (1) 

where a^(d) and h^{h) are the creation(annihilation) operators for vibrational phonons and cavity field photons, 
respectively. Eigen states of Pauli operator dz model the internal states of the j*'* two-level ion (j = 1,2,3) 
with energy splitting hwQ. We define the total spin operators as au = ^7=1 '^fc j where k = (z, -I-, — ) and use the 
eigenvectors of ct^ and d^ to represent the three ion internal states. The ionic internal states in product basis are 
labelled as ai Vi 'a), ' \ {ai' — —1 or +1)). The coupled basis vectors are |(t, cr^), where the label a{= 2s) refers to 

the eigenvalue of a^ given by aia + 2). The computational basis states, on the other hand, read as |ii, 12, 13) , where 
ij = for atom in ground state and ij = 1 for atom in the excited state. 

The interaction of cold ions, located at the node of the standing wave, with the quantized cavity field is given by 



Hi = hg{a+ + a-){V^ + h) sml'qia) + a)], 



(2) 



g being the ion-cavity coupling constant and 77 the Lamb-Dicke parameter. The interaction picture Haniiltonian 
obtained by applying the unitary transformation Uj = exp I —iHot/hj is a complex looking operator '30|. However, 

with the cavity coupled to red sideband of ionic vibrational motion {luq — uJc = v ), the relevant part of Hj in the 
rotating- wave approximation and Lamb-Dicke limit (77 ^ 1) reduces to 



Hii = hgri[(jjJba -\- aJw a)]. 

The possible values of a for three ions are 3, 1, 1, there being two distinct internal configurations that allow a 
The matrix T that transforms from the computational basis to the coupled basis is given by 
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The computational basis vectors are taken in the order |000), |100), |010), |110),|001), |101), |011), and |111), while 
the ordering of coupled basis vectors is |3, — 3), |3,— 1), |3,1), |3,3), |1,— l)i, |l,l)i, |1,— 1)2, |1,1)2- The subscript in 
|1, ±1)1^2 distinguishes the states with same value of a but different internal configurations. 

The product state of the composite system looks like \a,(Tz)®\'ni,n), where m is the number of center of mass 
vibrational quanta and n is the number of cavity field photons. The advantage of working in a coupled basis stems 
from the fact that the Hamiltonian of Eq. ^ does not connect states with different values of a. If the initial state 
is an eigen state of cr^ with eigenvalue ai{(Ji + 2), the system evolves into a linear combination of states of the type 
|tT7,crz,m,n) with —aj < a^ < Uj. 

For (T/= 3 the Hamiltonian (Eq. ^) in the basis |3, 3) (g) |m — 2, n — 2), |3, 1) (E)\m— l,n—l), |3, —1) (g) \m, n) , and 
|3, —3) IS) \m + l,n + 1) is written as 



iJ„(a = 3) 



/ V2A^n \ 

y/2Amn y^Bmn 
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(5) 



where 



and 



Defining 



Ar,in = ^gf]\\-^ijn- l)(n- 1),B,„„ == hgr]\j2mn, 



Cmn = ^g-n\j-z{m+ l)(n + 1). 



(6) 



(7) 



A* = (^m« + 5™„ + Cln), and /3 = Vm' - 4A^„C^, 
the eigenvalues of H ^ in units of hgrj , are 

El = -^fi-f3, E2 = Vf^-P, 



(8) 



The unitary matrix U that diagonahzcs H^^{a = 3), {Hjj{a = 3)|(^i) = Eilipi), (i — 1,4)) is given by 
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where A^i = M - 2A2„„ and ^2 = A* - 2C^„. 
Starting from a given initial state, 

is obtained by applying the transformation exp i—iHjjt/h]. An internal state of ions may be rewritten in terms of 
the coupled basis states as 



Cz iJz (Tz \ (g) \m, n), of the system, the time evolution in interaction picture 



r(l) ^(2) ^(3)\ 



,-i^\-f . = E4i-'-^ 



(10) 



The state |cr, CTz) iE" |to, n)^ of the composite system with m vibrational quanta and n photons may, in turn, be expanded 
in terms of eigen states of Hu (a) as 



E^i(^)i'^fc) 



(11) 



We recall that the matrix T operates on the internal states of the two-level atom, while U operates in the composite 
Hilbert space formed by internal states, cavity field states and the state of vibrational motion. The state of the system 
at instant t is given by 



*w = EE4^^(^)e^p(-^^fc5^^) i-^fc) 



(12) 



To go back to the computational basis we use the inverse transformation. The state obtained is in the interaction 

picture and we go to Schrodinger picture by the transformation ^s(0 = e B~^(i). As the state of the composite 
system evolves, the internal states of ions, the cavity field state and the state of vibrational motion become entangled. 



III. THE COMPOSITE SYSTEM STATE 



A. Initial state, *™+i,„+i(0) = |000) 0\m+l,n+l) 



Consider the three ions in their ground states with the center of mass motion cooled down to a state with m + 1 
vibrational quanta. Photon field set to a fixed detuning is injected in to the cavity such that the cavity state is an 
n + 1 photon Fock state at i = 0. The analytical expression for the state of the system at current time t is found to 
be 



*m+i,«+i(t) = aQ{t)\000)(g)\m+l,n+l)+ai{t)\Wi)(g)\m,n) 

+a2{t))\W2) ^\m-l,n-l)+ a3(t))|lll) ® |m - 2, n - 2) 



(13) 



where \Wi) = |ioo)+|ooi)+|oio) ^^^ |^^^ ^ |iio)+icn)+|oii) ^^^ ^ ^^^^^^ ^-^^^ ^^^ ^^^ ^^^ excited ions, respectively. 
The coefficients ai{t),i = to 3, 



ao(t) = {{13 - H2) cos {\/(p + f3)gr]tj + (/? + ^2) cos (^\/(p^~f3)gT]tjj , 



(14) 



ai(t) = -- 



^ — ILOlt 



2/3 



v/(/3-M2)(/3 + Aii)sin ^{fi + l3)gr^t 



+ V(/3 + M2)(/3-Aii) sin ( ^(p - (3)gf]t^ ) 



(15) 



a2(i) = ' '"'*^^' ^^ (cos [Vii:^W)g¥) - cos (^Vii^^^)9¥)) , (16) 



and 



asW = ^g- (V(/5 - M2)(/3 - Ml) sin (^^(m + 0)9Vt) 



2/3 
-\/(/3 + M2)(/9 + m)sin (v/(M-/3).9^i)) , (17) 

satisfy the normalization condition ^^^P |fli(i)l = 1- The frequency Wi = J^(?7i+|)+Ci;c(n+1) — -2|^, refers to zero point 



energy of the initial state. The probability amplitude a2(i) is zero whenever cos I \/(/i + (i)gTjt 1 — cos I \/(/x — (3)grjt 

For values of t such that cos I \/{p + i3)grit\ =cos ( \/{ii~ P)g7]t] = ±1, the composite system state is a separable 

state. When the condition sin ( ■\/(/x + I3)g'qt\ —sa\ ( •y/(/i — I3)grit\ = ±1 is satisfied, the composite system is found 
to be in the state 

*J^;i,„+i(t) = ai(t)ll^i) ® \m,n) + a3(i))|lll) ® |m - 2, n - 2), (18) 

where the probability amplitudes 

ai(t) = ^ ( V(/3 - M2)(/3 + Ml) + V(/3 + Ai2)(/3 - /ii)) , (19) 

az(t) = ^ (v/(/3 - M2)(/3 - Ml) - V(/3 + M2)(/3 + Mi)) , (20) 

depend strongly on the initial state photon and phonon number. The three ions in W\- like state are found to be 
entangled to photon-phonon state, constituting a state having four-partite entanglement. 

The probabilities Pi (t) of finding i number of ions (i = to 3) in excited state are plotted as a function of variable 
T in figure ([T]), for the choice to = n = 2. The value of parameter r = gr\t is determined by the cavity ion coupling 
strength, Lamb Dicke parameter and the interaction time. We notice that for r w 37r/4 the system is found to be in 
a separable state and for r « Stt/S in state Vl/g 3 , with P\{t) k, 0.75, and P-iij) « 0.25. 

B. Initial state, $m-2,n-2(0) =)|111> ® jm - 2,n - 2) 

When all three ions are in their excited states at i = 0, the center of mass prepared in a state with m — 2 vibrational 
quanta, and the cavity in n — 2 photon Fock state, the state of the composite quantum system at current time is 
found to be 

«'m-2,„-2(i) - a3{t))\000)(g,\m + l,n+l)+a2it))\Wi)(g,\m,n) (21) 

+ai{t))\W2) (» |m - l,n - 1) + ao(f))|lll) (» |m - 2,?i - 2). 

The maximum number of coupled basis states populated by the ions as the interaction time increases is four, 
independent of the initial state photon or phonon number. Recalling that for interaction time such that 

sin I \/{fi + I3)grjt\ = sin ( a/(m ~ 0)gvi) ^ ilj the probability amplitude ao{t) ~ a^^) — 0, we find the three ions in 

VF2-like state entangled to photon-phonon state such that 

^ZU,n-M = 03(0)1000) ® |m -H 1, n + 1) -f ai(i))|W^2) ® |to - 1, n - 1), (22) 

with probabihty amplitudes a\{t) and 03 (i) given by Eqs p9|) and ((20|) respectively. Figure ([2]) displays the probabilities 
Pi(j) (i = to 3) as a function of variable r for the choice m = 2,ri = 2 at t = 0. The initial state preparation in 
this case involves cooling the center of mass motion to zero phonon mode with ions in ground state and the cavity in 
vacuum state. The internal states of three ions can now be prepared in the excited state through controlled resonant 
interaction with external lasers. The ion's state is well preserved due to the long lifetime ('--^l s). We observe that the 
composite system periodically returns to initial separable state with a period of i « ^^ . The black arrows point out 

the Pi{t) corresponding to state ^^^(1). 



C. Quantum state control and the number of vibrational quanta 

The maximum number of coupled basis states populated by the ions can be controlled by the number of initial 
state photon and phonon number. With the center of mass prepared initially in one phonon mode (m +1 = 1), and 
the cavity having one or more photons (n + 1 > 1), the composite system is found to be in the state 

*i,n+i(i) = cos(V3(n+l)5^t)e-^'"^*)|000) ® |l,n+ 1) 

- i sin (V3(n+l)577t) e-''^'')\Wi) (g> |0, n), (23) 

whereas for ionic center of mass in zero phonon mode, the state of the system remains unchanged. The minimum 
interaction time needed to get the three ion Wi state generation probability peak is imin — / Deterministic 

2gr)^3(n+l) 

coupling of ion's quantized vibration in the trap to the cavity mode, has been demonstrated by Mundt et al. [l3|. 
For cavity ion coupling strength g = 8.95 MHz, Lamb Dicke parameter value rj = 0.01 and cavity prepared in single 
photon state at t = 0, we obtain imin — — ^^ = 10.133/i sec. This is the simplest setting for generating three ion 
W state with a single ion in excited state. The three ion \Wi) state generation time can be decreased by increasing 
the number of initial state photons in the cavity. The increase in cavity decay probability with increase in photon 
number is likely to reduce the \Wi) state generation probability and must be carefully accounted for. 
For the initial state preparation with m + I = 2 and ?i + 1 > 2 the state of the system at current time is 

^2.«+i(t) = ^ ((371 + 3) cos [v/2(5n + 3)ffr;^j + 2n j |000) |2,n + 1) 



.e-^&±ism 



(5n + 3) 



(v/2(5n + 3)577t)|M/i)®|l,n) 



+ ^^) (c°" (v/2(5n + ■S)mt) - Ij \W2) ® |0, n - 1). (24) 

The period for the system to return to initial separable state is t^ — , ^'^ At L„ = , ^ we have ions 

in W2-like state coupled to photon-phonon system in state 



^2,»+i(t) = '7.'*^r.+ ^^ |000)^|2,n+l) + ^''"7/!lr ^'^^)^'0'"-^ (25) 

(5n + 3) (5n + 3) 

We notice that the maximum probability of finding the state \y/-2) increases with n, approaching = || in the large 
n limit. When initial state phonon and photon number is greater or equal to three, all the four coupled basis vector 
states seen in cq. (IT^ can be reached. 

D. W- state Generation probabilities 

Three qubit W-states are extremely useful for implementing various communication protocols, quantum state 
transport and quantum gates. These states do not have genuine tripartite entanglement, but maximal pairwise 
entanglement. For these states, global negativity and partial two way negativity is non zero, while the partial three 
way negativity is zero [3l|, which is equivalent to the three tangle [32l | being zero. The tripartite entanglement of the 
state is due bipartite correlations. We label the three ions as subsystems A, B, C, and consider the phonons and 
photons to constitute a single quantum system D. The reduced state operator for the ions 

Pil^n+lit) = tro (|*™+l,„+l(t)) (*™+l,„+l(t)|) 

obtained by tracing over the vibrational and cavity state degrees of freedom is a mixed state given by 

Ptl'ln+iit) = l«o(t)r |000)(000| + |ai(t)|2 \W,){W,\ + \a^{t)\' \W^){W^\ + \a,{t)\' |111)(111|. (26) 

The probability Pi{t) and P2{t) in figure (jT]), are the probabilities of finding the three ions in the state \Wi) and 

37r 
8S»?' 



IW2) for the choice m + l = n + l = 3. We notice that for t = ^, Po{t) = P2 (''') = 0, Piij) shows a peak and Pz{t) 



is finite. 




FIG. 1: The probabilities Pi, i = to 3 versus r(= grjt) for the initial state |000, 3, 3). 



The choice m = n = Q, in Eq. p^ yields 

p^f^'it) = cos^ (VSgrit) |000)(000| + sin^ (^3^*) |VFi)(iyi|, 



with deterministic \Wi) state generation at t — /- , where k is an odd integer. 
For the choice to = n = 1, in Eq. (|13p we get 

P^^'^it) = (J cos {^g^i) + i^ |000)(000| + ^ sin2 (4^^) |W^i)(W^i| 
+ l(cos(4OTi)~l))2|M^2)(M^2|. 



(27) 



For an interaction time of i = ^^, fc = 1, 3, 5, ..., the ionic state \W\) is found with a probability of 75% in the state 



^ABC 



r^ui, 



p-i-r^ ^vSotJ " 16'°°°^^°°°' + Jiw^i)(w^il + ^|w^2)(M/2| 



(28) 



whereas the probability of finding the three ions in state \yV-i) is maximized for t = -^^ with the reduced state 
operator reading as 



fcTT 



P22 ( ^j = 4 1000) (0001 + -^W^){W^\ 



(29) 



For f = ^2L with fc = 0,2,4, ... the three ions are found in their ground states. The three ion state is a mixed state, 
i|Wi)(W^i| + ||W2)(W^2|, for g-qt values such that cos (Agrji) = -5. 

The three ion state operator (Pm^2'„_2(0 = ^^-D {\^m-2,n-2{t)) ($m-2,n-2(i)|)), reads as 



PX^%n-2{t) = |a3(i)r |000)(000| + \a,{tr \W,){W,\ + |ai(i)r \W,){W,\ + ^1^ |111>(111|. 



(30) 



For m = n = 2, the W2 state population probability at peak value is found to be 75% as seen in figure (2). The 
probability of populating Wi state is, relatively, small. 




270 



FIG. 2: The probabilities Pi, i = to 3 versus r(= grjt) for the initial state: |111, 0, 0). 



IV. ENTANGLEMENT OF IONIC QUBITS AND PHOTONS 



The entanglement of ionic memory qubits in W-like state to cavity photons is an extremely interesting and useful 
aspect of the proposed scheme. Entangled cavity photons can transport information to a remote cavity in a fast 
and reliable way. For analyzing the entanglement dynamics, we consider the (rn + 1 — i) phonons and (n + 1 — i) 
photons to constitute a single quantum system in the space spanned by vectors \i) (i = to 3), where \i) represents 
the state \m + l — i,n+l~ i). Furthermore, the ground and excited state of an ion represent logical bits |0) and |1) , 
respectively. The composite state of Eq. (fT3|) . in logical basis, reads as 



* 



m+1,^+1 



(t) = 



ao(i) 10000) +ai{t) 
+a2{t) 



|1001> + 10101) + 10011) 



x/3 
1102) + |1012) + 10112) 



a/3 



a3(t)|1113). 



(31) 



Labelling the three ions as A, B, C and cavity field plus phonon state as subsystem Z?, the possible bipartite partitions 
of the system are A - BCD, B - ACD, C - ABD, D - ABC, AB - CD, AC - BD, and AD - BC. We notice 
that the ionic state is symmetrical with respect to interchange of a pair of qubits. As such distinct partitions are 
reduced to A — BCD, D — ABC, and AB — CD. The state is in bi-orthogonal Schmidt form for qubits A, B, C 
and subsystem D. Quantifying or even detecting the entanglement of the composite system with four subsystems is a 
fairly complex task. The four subsystems may have four-partite GHZ state like correlations, or 4-partite entanglement 
resulting from bipartite entanglement between the subsystems. On tracing over anyone of the subsystems, 4-party 
GHZ like entanglement is destroyed, but the three remaining subsystems may have tripartite GHZ like entanglement 
or W-like entanglement. While the GHZ-like tripartite entanglement is completely destroyed on tracing over one 
subsystem out of the three, the reduced bipartite mixed state may still have bipartite entanglement (provided that 
the tripartite system had W-like entanglement). The entanglement distribution in the pure state of the composite 
system determines the entanglement available simultaneously to all the four parties, three selected parties or a pair 
of subsystems. In particular, the nature, distribution and dynamics of entanglement involving the memory qubits 
(trapped ions ABC) and photons is of special interest if the system is part of a quantum network. 

A. Peres [24l| was the first to observe that a partial transpose of density matrix associated with a separable state of 
a bipartite system is still a valid density matrix and thus positive (semi) definite. Horodecki et al. [25] proved that a 



positive partial transpose (PPT) was a necessary and sufficient condition for separability of a state if the dimension 
of the Hilbert space does not exceed six. A negative partial transpose is a clear signature of entanglement. Negativity 
|2ll . \22 . |23| based on Peres Horodecki PPT criterion has been shown to be an entanglement monotone [2a]. In a 
recent article |28l|. we introduced 2— way and 3— way negativities to discuss the entanglement of three qubit states. 
In refs. (27|, |29|, a characterization of multipartite quantum states having N subsystems, based on partial K—way 
negativities has been proposed. The /f— way partial transpose with respect to a subsystem is constructed from the state 
operator by imposing specific constraints on the matrix elements involving the states of K subsystems of multipartite 
composite system. It has been shown that the partial transpose of density matrix of an N-partite system, with respect 
to a given subsystem, can be written as a sum of i^— way partial transposes where 2 < K < N. Contribution to 
negativity due to a K—way partial transpose is easily calculated. The K—way negativity (2 < K < N), defined as 
the negativity of K—way partial transpose, quantifies the K—way coherences of the composite system. In contrast 
with the entropic measures of entanglement, where a reduced state operator is used to obtain information about 
the correlations between parts of a composite system, the partial K—way negativities are calculated from the state 
operator of the composite system itself. The negativity of global partial transpose remains invariant under local 
operations and classical communication (LOCC), however, the partial K—way negativities may increase or decrease 
under LOCC at the cost of each other. The partial K—way negativities show the entanglement distribution in different 
parts of the system in a given state. While the total entanglement in a composite system can not be increased by 
local operations, the entanglement distribution can be changed by local operations. In this section we investigate the 
entanglement dynamics of the composite system state by using the global negativity to detect the entanglement, and 
partial K—way negativities to determine the entanglement distribution amongst different subsystems. 

A. Definition of Global and K—v/ay negativities 

The global partial transpose p^* of four-party state operator p^^CD _ |\i/^ (\|/| with respect to subsystem A is 
constructed by using the condition 

{iii2i3i4\ Pg"^ \31j2j3j4) = (ji«2i3«4| p-^^"^^ \i1j2j3j4) ■ (32) 

The K—way partial transpose (2 < X < 4) of state operator p^^CD -^^j^j^ respect to subsystem A is defined as 

N 



{iii2i3U\ P^K \hhj3j4) = (jiJ2i3«4| p-^^"^^ \i1j2j3j4) if X! (-^ ^ ^ir^dm.) = K, 



(33) 



and 



771 — 1 



N 



(Z1Z2«3J4| Pk"" \jlj2j3j4) = («1«2«3«4| p"^^^^ \jlj2j3j4) if X! (^ ~ (^i™ J™) ^ ^> 



where 



^ir^Jm = 1 for «77l = J77 

^i^Jm = for im ^ j„ 



(34) 



Similar constraints are applied to construct global and K—way partial transpose with respect to subsystems B, C or 
D. The negativity of p^p is related to the trace norm of p'^p through 



NP 



dp- 1 



1 



(35) 



Negativity based on Peres Horodecki criterion is a natural entanglement measure. We define the global negativity 
with respect to subsystem p as 



"S = i 



hi 



Pg 



1 



(36) 



and K—way negativity [27|, |29| as 



""^dhd 



Pk 



(37) 
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where dp is the dimension of the Hilbert space associated with subsystem p. 

Furthermore, we define 3— way partial transpose p^^~^^'^ involving the qubits ABC as 

4 
if ^(l-(5,„j„J = 3,andii7^ji, (38) 

771—1 

{iii2i3H\p^''^'"^ \J1J2J3J4) = {Hi2i3H\p'^^^°\jiJ2J33i) 

4 

if ^i'^-Shr^Jr^) == 3,andii=ji, (39) 



and 

' - ^^"^^ IJ1J2J3J4/ = V«l'2*3«4|P' \J1J2J3J4:) n 

m— 1 

with the corresponding negativity 



»1*2*3*4| fK^^^^ IJ1J2J3J4) = (zii2Z3i4| P^'"''" IJ1J2J3J4) if ^ (^ " ^'™.^™) ^ 3, (40) 



j^A-ABC ^ - 
3 dp-1 



F3 



-1 . (41) 



Analogous definitions hold for the 3— way partial transposes p^^ abd ^ ^^a acd ^ ^^^ ^^d bco ^ 

B. Global and partial i("— way negativity 

Global negativity with respect to a subsystem can be written as a sum of partial X— way negativities. The global 
transpose with respect to subsystem p, written in its eigen basis is given by 

p^^=^Af+|vI.f+)(*f+|+^AM*^)(^f-|, (42) 

i i 

where A^ ~^and |^^ ^) (A,- "and |^,- ~)) are the positive (negative) eigenvalues and eigenvectors of pj", respectively. 

Pq 1=1, the negativity of p^ is given by 

^s=-^E(*M?ai*n=-^E^r- (43) 

The global transpose with respect to subsystem p, may also be rewritten as 

N 

fj = Y.fK-(^-^)P- (44) 

K=2 



Substituting Eq. dS]) in Eq. (gS]), we get 

N 

-2EAf-=-2EE(*M^1*f")+2(^-2)E(*M^I*f")- (45) 



i A'=2 i 

Defining partial iiT— way negativity E^ {K = 2 to N) as 



and 



^]^ = -^E(*M/ti*r>' (46) 



^? = -^^^E(*M?l*f"> (47) 
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Initial State |0000> 




FIG. 3: The global negativity Nq , and entanglement measures Ej^, for K = 2 to 4 versus t{— g-qt) for the initial state |0000). 



we may split the global negativity for qubit p as 



N 



nI = Y.Ek~ei 



K=2 



(48) 



To obtain tripartite GHZ state like correlations between three subsystems, we calculate 

2 



E. 



A-ABC 



tE(* 



r3 \ i 



(49) 



E. 



A-ABD 



TE(*M?^'"'"°i*r)' 



(50) 



and 



TT^A-ACD 
^3 



dp- 



tE(*M?3 ^"1*^)' 



(51) 



It is easily verified that for the system at hand £"3 
qubits have no genuine tripartite entanglement. 



A-ABC 



= 0, and E. 



A-ACD _ pA-ABD 



E^/2.0. As such the three 



C. Entanglement dynamics of the pure state "f^ 



L(i) 



We use the global negativity to detect the entanglement of parts in bipartite splits of the system. In case the 
negativities Nq, Nq , Nq , Nq and Nq^ are non zero, the system has 4-partite entanglement. The 4— way, 3— way 
and 2— way partial negativities identify and quantify the different types of entanglement between the subsystems 
of the composite system in a given state. From the state operator p{t) — \'im+i,n+i{t)) (^m+i,„+i(t)| , the global 
and ii'— way partial transposes are constructed by following the prescription given in section IV ( Eq. (|32p and Eq. 
([55)) ). The state of Eq. P^ is very special in that analytical results can be obtained for partial K—way negativities 
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Initial State |0000> 




FIG. 4: The global negativity Nq , and entanglement measures E'j^, for if = 2 to 4 versus r(= gr}t) for the initial state |0000). 

characterizing the state. Using negative eigen functions of p^ ^ the entanglement measures i<^^ for K ~ 2^i and 4 are 
easily obtained. The negativity of global partial transpose with respect to qubit A is found to be 



N^ = 2 



\ 



laowr 



2|«iW|' , \a2(tf 



|a3(t)| 



2 , |ai(t)|' , 2\a^{t)f 



Using the eigenvector corresponding to the negative eigenvalue of p^^ , we get the partial negativities 



(52) 



Et = —\\ao{tr\a,{tr + 



2 , laiWl'NWr 



(53) 



j^A^J_ 2\ao{tf\a2{t)\' 2\a,{t)\'\a,it)\' 
^ Ng\ 3 3 



(54) 



and 



E^^ 



Ng 



\ao{t)\'\a,{tt , \a2it)\'\a,it)\' ^ J \aiit)\' ^ \a2it)\' 



To J'd J^d 



„Td 



(55) 



It is easily seen that Nq — Nq = Nq . Next we construct the transposes Pq" , p^" , ^3° , and P2° , for the decomposition 
ABC — D oi the state ^m+i,ri+i(0 ^^'^ obtain 



NS = -\a,{t)\{\a,{t)\ + \a2{t)\ + \a,m 

+ I\a^it)\i\a2it)\ + \a,m + 2\a2it)\\a,it)\ 



(56) 



E. 



D 



\Mt)\\a3{t)\ + -\a,it)\\a2it)\, 



(57) 
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FIG. 5: The global negativity Nq , and entanglement measures Ej^, for K = 2 to 4 versus t{— g-qt) for the initial state |1113). 



E^^^\ao{t)\\a2{t)\ + -\a,{t)\\a:,{t)\, 



(58) 



and 



E? = 3 l«oWI \ai{t)\ + - \a2it)\ |a3(t)| + - |ai(t)| \a2{t)\ 



(59) 



Since the three qubits have no genuine tripartite entanglement, the partial negativity E^ represents the entanglement 
of quantum system D with the W like states of three qubits. 

Treating AB as a single system in Hilbert space of dimension four, the global negativity is found to be 



Nq^ = - (/iO/Ui + /iO/U2 + /ilM2) , 



(60) 



where 



Mo 



i«o(i)r 



i«iwr 



Ml 



'2|ai(i)|' , 2|a2(i)|' 



M2 



i«3(i)r 



i«2(i)r 



(61) 



Figures (3) and (4) display the global negativity and K—way entanglement of qubits A and D in the state \E'33(i), 
as a function of interaction parameter r = g?]t. At r « ^, we have E^ = E^ = 0.5, indicating that the qubit A (or 
-B or C), has equally strong bipartite and tripartite correlations. Recalling that no genuine tripartite entanglement 
exists between the three qubits as evidenced by E^~ = 0, we have here the three qubits in Wi-like state entangled 

to the subsystem D. In other words, the cavity field is entangled to three ions in VFi-like state and can be used to 
transfer the entanglement of the composite system to a remote quantum system. The probability plot of figure ([1]) 
confirms that for interaction time t = -^^ -fo(''') = P2{t) = 0, Pi{t) shows a peak and Psir) is finite. A measurement 
that finds the cavity in two photon state collapses the composite system state to three ions in VFi -state with center 
of mass in two phonon state. 
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FIG. 6: The global negativity N§ , and entanglement measures E^, for K = 2 to 4 versus r(= gT]t) for the initial state |1113). 



D. Entanglement dynamics of the pure state $ 



m — 2,n — 2 



it) 



Analytical expressions for negativities and partial K—way negativities for the state ^m-2,n-2{t) (Eq. (|2ip ). written 
in computational basis as 



$ 



m-2 



,„_2W = a3{t)\0000) + a2{t) 



/ |1001) + 10101) + 10011) 



V3 



,^,^ 11102) + 11012) + 10112) ^ , 
+ai(i) I ^= I +ao(i)|1113). 



(62) 



are analogous to those for the state *m+i,n+i(i)- For the special case of state $oo(i), the global and partial K—way 
negativities are displayed, in Figs. (5) and (6) for qubits A and D, respectively. The negativities Nq^ for the states 
\I'33(i) and $00 (i) are shown in Fig. (7). At r = ^, the coefficient 02 (i) = ao{t) = , and E^ = E^ ^ 0.5 giving 
(using Eq. [22|) the state 

Svr 
*oo(t-== — ) = a3(i)|000,m + l,n+l) + ai{t)\W2,m- l,n- I), 



In this case we have three ions in W2-like state entangled to photon-phonon system. As such detecting the cavity in 
single photon state ensures that the three ions are in W2 state. The entanglement of the cavity field with three ions 
in VF2-like state may, on the other hand, be used to communicate with a remote quantum system. For an interaction 
time t = ^, P2{t) shows a peak and Po{t) is finite in the the probability plot of figure ^. 

The number of initial state vibrational quanta also controls the nature of entanglement in the composite system 
states. We notice that for a single phonon initial state the general form of the state (Eq. (P5|) ) 



*i,n+i(i)=ao(i)|0000)+ai(t) 



1001) + 10101) + 10011) 
V3 



(63) 



allows only bipartite entanglement. For m > 2, the composite system may have genuine 4-partite, tripartite as well 
as bipartite entanglement. This is an interesting aspect unique to systems where vibrational motion of ions is coupled 
to cavity field. 
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FIG. 7: The global negativity Nq^ versus r(= g-qt) for the initial states |0000) and |1113). 



WHAT DO THE 7^ WAY NEGATIVITIES ADD TO OUR KNOWLEDGE OF THE SYSTEM? 



Figures (3-7) show the entanglement distribution between possible entanglement modes available to the system in 
states ^33(7) and $00 (i)- Since analytical expressions for partial K—wsy negativities are available for specific initial 
state preparations, the rate of change of partial K~w&y negativity for a particular mode can be easily obtained, if 
needed. The two way partial negativity is seen to grow at the fastest rate being the first to reach it's peak value, 
followed by two peaks showing maxima of partial 3— way and 4— way negativities, respectively. Besides that a reversible 
entanglement exchange between different entanglement modes is observed. In Fig. (3) at the maxima of partial 4— way 
negativity, partial 2— way and 3— way negativities are rather small, whereas the minima of E^ correspond to a large 
contribution to total entanglement from bipartite and genuine tripartite entanglement. Similar trend is seen in Figures 
(4-6). Figure (7) complements the information about entanglement distribution obtained from Figures (1-6). Once 
the interaction is switched on, a typical qubit pair is found to be in an entangled state until the composite becomes 
separable again. As seen from the global negativity plots, the period after which the composite system becomes 
separable is given by T = j^. 

In particular, we look at the entanglement of states ^jy^\in+i (Eq- CHI)) and '^rn-2 n-2 (^q- (|22ll ). For both types 
of states A^^ = E^^ , that is there no genuine four partite entanglement amongst the subsystems A, B, C and D. The 
global negativities are however finite for A, B, C, D as well as AB, pointing to four partite entanglement. Although 
for qubits A, B, and C we get finite partial 2— way negativities, four partite entanglement cannot be due to two-way 
correlations because i?f = 0. The 2— way negativity of the partial transpose of Wi state or W2 state is 0.94, where 
as the three way negativity is zero. As such a W state is a state with maximal tripartite entanglement generated 
by 2— way correlations. In states ^„+i „_|_i and ^Z.-2 n-2 i^^s in a ly like state are entangled to subsystem D. We 
conclude that the four-partite entanglement of three ions with photon-phonon system is generated by 2— way and 
3— way correlations. 



VI. 



CONCLUSIONS 



We have studied the entanglement dynamics of spatially separated three two-level cold trapped ions in a high 
finesse cavity with the cavity tuned to the red sideband of ionic vibrational motion. Analytical expressions for the 
state of composite system as a function of interaction time are obtained for interaction Hamiltonian of Eq. [3] with 
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the cavity prepared initially in different photon number states. With three ions initially in their ground state, the 
state |000,77i + l,n + 1) evolves (i) for m = and n > 0, in a bi-dimensional subspace , (ii) for m = 1 and n > 1 
in a tridimensional subspace, and (iii) for m > 2 and n > 2 in a four-dimensional subspace of the coupled basis 
ionic states. The initial state |lll,m-|- l,n-|- 1) always evolves in a four-dimensional subspace independent of the 
initial phonon and photon number. The number of initial state vibrational quanta offers a control mechanism for 
manipulation of composite system states, when ions are in their ground state initially. This is an interesting aspect 
unique to systems where vibrational motion of ions is coupled to cavity field in contrast to the ions coupled only to 
quantized cavity field [33|, [34] or only with the vibrational modes 35]. Useful practical applications to implement 
information processing and communication related tasks, can benefit from this special feature. The reduced three 
ion state operator is obtained by tracing out the phonon and photon degrees of freedom. For cavity ion coupling 
strength g = 8.95 MHz, Lamb Dicke parameter value rj = 0.01 and cavity prepared in single photon state at i = 0, the 
minimum interaction time needed generate a three ion W-state is found to be ~ 10.133^ sec. For the initial state in 
which three ions are in their ground state and center of mass in two phonon state, the W2 state generation probability 
increases while the minimum interaction time to get the probability peak decreases, with increase in the number of 
photons present in the cavity at i = 0. 

The ionic qubits in W state are found to be entangled to cavity photons, that may be used to transport information 
to a remote cavity in a fast and reliable way. Multipartite entanglement dynamics of the composite system is examined 
using global, four-way, three-way and two-way negativities. Analytical expressions for partial K—way negativities 
{K = 2 to 4) are obtained. For the three ions prepared initially in their ground state or in their excited state, the partial 
K—way negativities are calculated numerically and plotted as a function of interaction parameter. These plots show 
the entanglement distribution as well as the rate of change with time of entanglement between possible entanglement 
modes available to the system. Besides that reversible entanglement exchange between different entanglement modes 
is observed. For specific values of interaction parameter, the three ions and photon-phonon system are found to 
have four partite entanglement, generated by 2— way and 3— way correlations. Three ions in W state are found to 
be entangled to photons. We expect this analysis to add to the understanding of multipartite entanglement in the 
context of trapped ions interacting with photons in optical cavities. 

Acknowledgments 

S. S. Sharma acknowledges financial support from FAEP/UEL, and Fundagao Araucaria PR, Brazil. E. de Almeida 
thanks Capes, Brazil for financial support. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 

[7 

[8 

[lo; 

[11 
[12 

[13 

[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22' 
[23 



A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 

C. H. Bennett and S. J. Wiesner, Phys. Rev. Lett. 69, 2881 (1992). 

C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres, and W. K. Wootters, Phys. Rev. Lett. 70, 1895 (1993). 

C. Monroe, D. M. Meekhof, B. E. King, W. M. Itano and D. J. Wineland, Phys. Rev. Lett. 75, 4714 (1995). 

D. J. Wineland, C. Monroe, W. M. Itano, D. Leibfried, B. E. King, and D. M. Meekhof, NIST J. Res. 103, 259 (1998). 
C. F. Roos, Th. Zeiger, H. Rohde, H. C. Nagerl, J. Eschner, D. Leibfried, F. Schmidt-Kaler, R. Blatt, Phys.Rev. Lett. 83, 
4713 (1999). C. Monroe, D. M. Meekhof, B. E. King, W. M. Itano and D. J. Wineland, Phys. Rev. Lett. 75, 4714 (1995). 
H. Haffner, W. Hansell, C. F. Roos, J. Benhelm, D. Chek-al-kar, M. Chwalla, T. Korber, U. D. Rapol, M. Riebe, P. O. 
Schmidt, C. Becher, O. Giihne, W. Diir and R. Blatt, Nature 438, 643 (2005). 

J.I. Cirac et al., Phys. Rev. Lett. 78, 3221 (1997). 

V. Buzek, G. Drobny, M. S. Kim, G. Adam, and P. L. Knight, Phys. Rev. A56, 2352 (1997). 

S. Shelly Sharma, Phys. Lett. A 311, 187 (2003). 

S. Shelly Sharma, N. K. Sharma and E. de Almeida, J. Phys. B: At. Mol. Opt. phys. 39 695(2006). 

A. B. Mundt, A. Kreuter, C. Becher, D. Leibfried, J. Eschner, F. Schmidt-Kaler, R. Blatt , Phys. Rev. Lett. 89, 103001 

(2002) 

A. B. Mundt, A. Kreuter, C. Russo, C. Becher, D. Leibfried, J. Eschner, F. Schmidt-Kaler, R. Blatt, Appl. Phys. B 76, 

117 (2003). 

J. Eschner, Ch. Raab, A. Mundt, A. Kreuter, C. Becher, F. Schmidt-Kaler, and R. Blatt, Fortschr. Phys. 51, 359 (2003). 

G. Morigi, J. Eschner, and C. H. Keitel, Phys. Rev. Lett. 85, 4458 (2000). 

C. F. Roos, D. Leibfried, A. B. Mundt, F. Schmidt-Kaler, J. Eschner, R. Blatt, Phys. Rev. Lett. 85, 5547 (2000). 

T. E. Tessier, I. H. Deutsch, and A. Delgado, Physical Review A 68, 062316 (2003). 

Tavis and F.W. Cummings, Phys. Rev. 170, 379 (1968). 

C. Di Fidio, S. Maniscaico, and W. Vogel, Phys. Rev. A, 65, 033825 (2002). 

D. Bouwmeester, A. Ekert, and A. Zeilinger, The Physics of Quantum Information (Berlin: Springer) (2000). 
K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewenstein, Phys. Rev. A 58, 883 (1998). 

G. Vidal, J. Mod. Opt. 47, 355 (2000). 

J. Eisert, PhD thesis (University of Potsdam, February 2001). 



17 

[24] A. Peres Phys. Rev. Lett. 77, 1413(1996). 

[25] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1(1996). 

[26] G. Vidal and R. F. Werner, Phys. Rev. Vol. 65, 032314 (2002). 

[27] S. S. Sharma, and N. K. Sharma, quant-ph/0608062 (unpublished). 

[28] S. S. Sharma, and N. K. Sharma, Phys. Rev. A 76, 012326 (2007). 

[29] S. S. Sharma, and N. K. Sharma, Phys. Rev. A 77, 042117 (2008). 

[30] S. S. Sharma, A. Vidiella-Barranco, Phys. Lett. A 309, 345 (2003). 

[31] S. S. Sharma, and N. K. Sharma, arXiv:0806.0887vl [quant-ph]. 

[32] V. Coffman, J. Kundu, and W. K. Wootters, Phys. Rev. A 61, 052306 (2000). 

[33] K. Fujii, K. Higashida, R. Kato, T. Suzuki, and Y. Wada, Int. Jour, of Geometric Methods in Modern Physics 1, 6 (2004). 

[34] JF. Cai and HP. Liu, Commun. Theor. Phys. 43, 3, 427 - 431 (2005). 

[35] A. Retzker, E. Solano and B. Reznik, Phys. Rev. A 75, 022312 (2007). 



